Implementation of many-qubit Grover search by cavity QED 
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We explore the possibility of A^-qubit (A'' > 3) Grover search in cavity QED, based on a fast 
operation of A''-qubit controlled phase-flip with atoms in resonance with the cavity mode. We 
demonstrate both analytically and numerically that, our scheme could be achieved efficiently to find 
a marked state with high fidelity and high success probability. As cavity decay is involved in our 
quantum trajectory treatment, we could analytically understand the implementation of a Grover 
search subject to dissipation, which would be very helpful for relevant experiments. 
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' Cavity quantum electrodynamics (QED) has got much progress over past several years People have achieved 
various optical and microwave cavities, with which some remarkable quantum phenomena have been demonstrated, 
' such as strong coupling of a single atom with a single photon f2|, quantum logic gating between two atomic qubits 
Hi photonic blocking [4], efhcient generation of single photons 5], and so on. 
In the present work, we will investigate how to realize a Grover search with more than three qubits by cavity 
^ QED. The Grover search, as one of the most frequently mentioned quantum algorithms, could effectively exemplifies 
^ , the potential speed-up offered by quantum computers 0. Recently, many authors have dedicated themselves to the 
qh" search algorithm by adiabatic evolution methods or nonadiabatic scheme \Td\. Although achievement of these 

I— I, schemes needs stringent conditions and demanding techniqu es they are really wonderful ideas. On the other 
hand, some authors had addressed the effect of decoherence [l2j|. noise [l3l |. and gate imperfection or errors [3] on 
[ the efficiency of quantum algorithms. We also noticed that, there had been intensive interests in achieving Grover 

• search al gori thm theoretically and ex p erimentall y by using NMR systems ^1^], trapped ions [16., 17, 18], linear optical 
— . [ elements |19l. l20j . cavity QED [U [l^, IH, [H, , and superconducting meso-circuits [11]. 

. Our interest is in the experimental feasibility of Grover search by a microwave cavity. A recent work [2^, with 
' two of the authors joined in, showed a three-qubit Grover search implementation in a microwave cavity. As it makes 
, use of the resonant interaction of the atoms with the cavity mode and consideres the detrimental influence from the 

• cavity decay, the implementation of [1^ is close to the reach of current technique. However, that scheme is hard to 
OO ! be extended to the case with more than three qubits due to the method [13] used there. In the present work, we will 

• show that, after some modification of the method in [13], we could accomplish a high-fidelity conditional phase flip 
^ [ (GPP) for many atomic qubits flying through a microwave cavity, based on which an A^-qubit (A^ > 3) Grover search 

might be achieved in a straightforward way. The main idea is that, we make a smart qubit encoding on the atomic 
levels, and send the atoms through a microwave cavity for resonant interaction. As long as we could exactly control 
the interaction time and the trajectories of the flying atoms, we may accomplish the A^-qubit GPP gate by one step of 
implementation, based on which we could efficiently carry out an A^-qubit Grover search with high fidelity. However, 
as the Grover search with more than two qubits works only probabilistically and also as our one-step implementation 
of many-qubit CPF gate intrinsically owns imperfection, the situation in our case is much more complicated than in 
two-qubit case. To make our scheme more realistic, we involve the cavity decay in our treatment, and as done in [23| . 
we will try to present some analytical expressions for the state evolution during our implementation. As a result, we 
could explicitly assess how well a Grover search is made subject to these disadvantageous factors. 
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2. N-QUBIT CPF GATE SUBJECT TO CAVITY DECAY 

As described in [2^ [23| , we consider the resonant interaction of N three-level atoms with a single- mode cavity. As 
shown in Fig. 1, the three levels can be described by \ij), \gj), and \ej), in which the subscript j means the jth atom, 
the states jcj) are excited states and \gj) are ground states. As \gj) and \ej) are coupled in resonance by the cavity 
mode, the states \ij) are not involved in the interaction with the cavity mode throughout our treatment. We encode 
the first qubit |0)^ (|l)i) on |ei) (|gi)), and other qubits |0)j. on jz^) {\gk)) with k = 2,3, ...,N. Considering 

the cavity with weak decay, which implies no photon actually leaking out of the cavity during our implementation, 
we may write the Hamiltonian in units oi h — I as. 



N 



(1) 



i=l 



where fli is the coupling constant of the ith atom to the cavity mode, and the atomic spin operators for raising and 
lowering are af' = \ei) {gi\ and = \gi) {ei\ , respectively, (a) is the creation (annihilation) operator of the cavity 
mode. We will first consider an implementation of iV-qubit CPF under cavity dissipation, based on which an 7V-qubit 
Grover search would be achieved in the next section. 

Assuming the initial atomic state in |ei) Y[i j=2 i=ij Idi) Ni) ^^'^ cavity initially in the vacuum state |0), we could 
deduce the evolved state by directly solving the Schrodinger equation regarding Eq. (1), which yields [27l |. 
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where G,,,. = ^] and Am+i,K = \lllT=i ^? ~ with m the number of the atoms initially in the ground 



state \g). 

In order to accomplish an A^-qubit CPF gate, we should also pay attention to the evolution with the initial states 



ei)n^2l«j) which is 
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with = a/IIJ — K^/IG. Assuming the coupling constant satisfies the condition = 1^2 = fia = • • • = rjjv ;:g> fii, 
we have, under the approximation of short time evolution with sin {Ai,i^t) << cos {Ai,K,t) , Eq. (3) reduced to 
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with a = 
to 



exp (— Kt/4) (cos {AiK,t) + sin (Ai^t) ) . Under the same condition as above, Eq. (2) similarly reduces 
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. In order to make a CPF gating, 
fll — k'^/16. As a result, the coefficient becomes 
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where i9 = 7r-\/l + IGmTy^/ (16 — y?), j] = Q/fli and a = — exp (j-nii/ ^16 — /i^^ with /i = K/fJi. It is eas- 
ily verified that a ^ — 1 and /3m — > 1 in the case of k 0. We may take /V = 4 as an example to show 
the CPF more explicitly. Following the deduction above, we have the four-qubit CPF written as U^^^ = Jqooo = 
dmg{a,/3i, /3i,/3i, /32,/32,/32, /33,1,1,1,1,1,1,1,1} with respect to the bases \e1i2i3ii), |eii2»354), {eiiigsii) , \e1g2i3u), 

\eii29394), \eig2i394), |eiff253«4), 161329354), |ffl«2«3«4), \9ii2i394), \gii293H), \9i92i3i4), 151*23354), |5l52«354), 

|5i5253«4), and I51525354) , respectively, where the latter eight states do not evolve under the Hamiltonian of 
Eq. (1) in the case of vacuum cavity state. To assess the CPF gating, we assume the atoms to be initially in 
l^'o) = (I51) + |ei)) (132) + 1*2)) • • • (|.9jv) + Vn))- We may calculate the fidelity and the success probability by 



employing the relations F = {-^ f\ {/("'t j^g) [/(») |^^) ps] and P = (*/|^'/) [23,123 As shown in Fig. 2, we could 
reach considerably high fidelity F and high probability of success P as long as 77 — > 10, even in the case of relatively 
large decay rate. In what follows, we will employ our CPF to carry out Grover search under weak dissipation. 



3. N-QUBIT GROVER SEARCH 



Generally speaking, provided that initial state of the system has been prepared in an average superposition state 



jvE'o) = (1/V2^) ^ \i) , Grover search can be depicted as the iterative operation D'^^^J^ by at least 7r\/2^/4 times for 

finding the marked state |t) with an optimal probability, where the quantum phase gate Jr = 1—2 \t) (t| (with / being 
the identity matrix) plays an important role to invert the amplitude of the marked state, and the diffusion transform 



D(^) is defined as = 2/ A' 
I51), N2), I52) 



5ij — 1,2,..., AT, AT = 2 ). Considering the qubit subspacc spanned by {|ei), 
«Ar_i), I5JV-1), I^Af), 15^)}, to make a Grover search, we could construct following transformation, 

Q{N) ^ w^^4^\w®^Jp = W^^/f^...^ W^^Jp = -D^^^Jp, (7) 



where JqqIq — diag{—l, 1, • • • , 1} = /^^^ — 2 |00 • • • 0) (00 • • • 0| (with superscript (N) denoting A^-qubit case and • • • 
denoting A" elements) and the Hadamard gate in our computational subspace is given by 
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In the case of A^ qubits, Eq. (7) implies that the diffusion transform Z)'^^ = ~T^'^^'^eii2 - i„^'^^ ^® always unchanged 
no matter which state is to be searched. The only change is the operator Jp for different marked states. Based on 

the CPF gate constructed in last section, we have J( 

the gate Jiili as. 



(AT) 
00---0 



diag{ — l, 1, • ■ ■ , 1} = ujj^p, from which we could achieve 
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where 



l«j)(5j 



\gj){ij\ with j 7^ 1, and a^^i = \ei){gi\ + |gi)(ei|. To achieve Q^^-*, we will con- 
struct the CPF gate Jp = I — 2\p) {p\ . For example, in the case oi N = 10, the number of the involved 
quantum states is 2^*^, and the operation is to label a marked state by Jp with p one of the states from 
{|0000000000),|0000000001),|0000000010),--- ,11111111111)}. To carry out the ten-qubit Grover search, we need 
two ten-qubit Hadamard gates W'^^'^. Based on the CPF gate marking the state |5i52535455565758595io) , we could 
construct other 2^" — 1 gates for the marking job by single-qubit rotations. For example, 



313233349536373839310 ~ "^1111111111 — 

Sx,10Sx,9Sx,8Sx,7Sx,6Sx,5Sx,4:Sx,3Sx,2O'x,lJ0000000000(^x,lSx,2Sx,3Sx,4:Sx,5Sx,6Sx,7Sx,8Sx,gSx.l0- 

(10) 
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So with a state marked, and with the ten-qubit diffusion transform l)^'^^) which is generated by combining two 
Hadamard gates with Jooooooooooj a full Grover search for ten qubits is available. 

Let us consider a ten-qubit Grover search for the marked state |eiZ25'35'4i5i6^7*8*9*io) as an example. The Rydberg 
atoms from the 1st to the 10th are input, as shown in Fig. 1, with principal quantum numbers 49, 50 and 51 denoting 
\gj), \ij) and \ej), respectively. Suppose the atoms initially prepared in a superposition j^Po) — x (|ei) + |gi))(|i2) + 
Iff2))(|«3) + I53)) • • • (|«io) + Isio))- We plot in Fig. 3 for the success rate of the Grover search in different cases. As 
the Grover search involving more than two qubits is intrinsically probabilistic, how to obtain an optimal search is a 
problem of interest. In contrast to the previous study for an ideal Grover search [29], our analytical expressions along 
with some numerics could identify the optimal search under cavity decay and other imperfection due to gating. We 
know from Fig. 3 that, because of the cavity decay, the more iteration steps, the more detrimental effect from the 
cavity decay involved. As a result, we prefer to have less iteration steps for an optimal search. The detrimental effect 
from the cavity decay is also reflected in the figure that the maximum value of P is slightly lower in the case of larger 
decay rate. 

4. DISCUSSION 

Our CPF gating plays the essential role in our Grover search scheme, which much reduces the implementation time 
compared to that by a series of two-qubit conditional gates and single-qubit gates [13] . Sec II has shown that the CPF 
gating time only depends on the weakest coupling Qi and the cavity decay rate k, but irrelevant to the qubit number. 
So our scheme could keep constant implementation time no matter how many atoms are involved, which favors a 
large-scale QIP. As a result, if we suppose the single-qubit rotation takes negligible time compared to the CPF gating, 
then our Grover search could be carried out by a constant time irrelevant to the qubit number. Specifically, assuming 

= 2-77 X 49 kHz [1], Qi = 2iT X 4.9 kHz, and k = O.lfli, we have to = ir/y/^l — « 102 us, which is much 

shorter than either the cavity decay time, i.e., 27r/K « 2 ms, or the Rydberg atomic lifetime 30 ms [l|. However, for a 
Grover search, the operation Q*^^^ has to be made for several times depending on the specific search task. In the case 
of iV = 10, the total required time for a Grover search would be about 2[7r\/2^/4] x io ~ 5 ms which is comparable 
to the cavity decay time. 

As shown in Fig. 3, our numerics shows that the iteration steps for optimal search are 25, 23, 17, respectively, in 
the case of = 0, 0.05, and 0.1. Since the operation time to = nf-^fll — k^/16 varies very little in the case of weak 
cavity dissipation with respect to the ideal case, the less iteration steps under dissipation would definitely lead to 
reduction of our implementation time. Straightforward calculation shows that the required time for an optimal search 
are 5 ms, 4.6 ms, and 3.4 ms, respectively, in the case of /i = 0, 0.05, and 0.1. It is of importance to have a fast 
implementation of Grover search in view of decoherence. But if we seriously consider the values in above calculation, 
we could find that fi = 0.05, and 0.1 correspond to dissipation time of 4 ms and 2 ms, respectively. As a result, to 
have a fast implementation with our scheme, we must increase the coupling strengths SI and f2i by at least ten times, 
which makes sure the required implementation time to be much shorter than the dissipation time. 

In current microwave cavity experiments [l|, we may consider the interaction between atoms and the cavity mode 
by ricos(27rz/Ao) exp(— r^/w^) ^ cos(27rz/Ao) [H, where VL is the couphng strength at the cavity center, r is the 
distance of the atom away from the cavity center, and Aq and w are the wavelength and the waist of the cavity mode, 
respectively. To meet the condition of = lOSli, we may send the first atom going through the cavity along y axis 
deviating from the nodes by (Ao/27r) arccos(l/10), but other atoms through the antinodes. Specifically, for the N 
input atoms (suppose N being an even number for simplicity), we should have the tracks of the atomic movement as 
zi = -{N/2)\o + (Ao/27r) arccos(l/10), Z2 = -[N/2 ~ l)Ao, Zn = {N/2)\o. Moreover, with current cavity QED 
techniques, the design in Fig. 1 could be realized by 50 separate microwave cavities with each Ramsey zone located 
by a cavity. Since each microwave cavity is employed to carry out a ten-qubit phase gate Jqq^Io^ cavities should 
be identical. While searching for different states, we employ different single-qubit operations. So the Ramsey zones 
should be long enough to finish the necessary single-qubit operations. 

However, it is highly challenging with current experimental technology to simultaneously send many Rydberg atoms 
through a cavity with precise tracks and velocities although two Rydberg atoms going through a microwave cavity 
simultaneously have been achieved [2l[. Here we assess the infidelity due to the imperfection in atomic velocity by 
a simple example below. Consider the simple situation that the first atom moves a little bit slower in the cavity 
than other atoms, that is, ti = to + St, t2 ^ t^ ~ ■ ■ ■ ~ i„ = to, with to = 7r/-\/g^ — k^/IG. As shown in Fig. 4 
with 77 = 10, we can see that the infidelity increases with both St and /i. In actual experiment, the situation is much 
more complicated. Besides the diversity of atomic velocity, there will be other imperfections such as classical pulse 
imperfection, differences of the cavities, deflected atomic trajectories and so on. Our simple estimate could show the 
importance of highly controllable movement of the atoms in implementing our scheme. 
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5. CONCLUSION AND ACKNOWLEDGEMENTS 

In conclusion, we have studied many-qubit CPF gate in cavity QED under weak cavity dissipation, and explored 
the possibility of many-qubit Grover search in a dissipative cavity QED system. As a theoretical work, like previous 
publication in cavity QED, e.g., [3l| . we have simplified the complexity in the realistic system, and tried to present 
some analytical expressions for the state evolution subject to the cavity decay. Although we took an example involving 
ten qubits, our approach could also be easily used for three or four atomic qubits in cavity QED, which might be closer 
to experimental realization. On the other hand, we have noticed experimentally achieved entanglement of six photons 
[3^ , eight trapped ions and twelve NMR ensemble qubits . So we also hope our idea could be extended to other 
systems, such as to trapped ions. Furthermore, from above discussion, we have known the experimental difficulty for 
ten atomic qubits to be entangled and used for Grover search by current cavity QED technique. Nevertheless, we 
argue that our study would be useful for future QIP experiments by microwave cavity. 

One of the authors (H.Q.F) would like to thank Professor Kelin Gao, Zhenyu Xu and Xiaodong He for helpful 
discussions. This work is partly supported by National Natural Science Foundation of China under Grants No. 
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Captions of Figures: 

Fig. 1. Schematic setup for finding the marked state |eii2.93.94i5i6^7^8*9Ho) in ^ ten-qubit Grover search, where ten 
atoms are simultaneously sent through the cavity with proper speed. The ten state-selective field ionization detectors 
£>i, • • •, Dio are settled at the end of the passage for checking the states of the qubits. The operations U^^^\ H^^^, 
Sx.3, and SxA are defined in the text. The inset shows the atomic level configuration, with bold lines for the states 
encoding qubits and the arrows for the coupling between the atoms and the standing wave of the cavity. 

Fig. 2. The variations of the fidelity and probability of success with respect to both rj and fj,. 

Fig. 3. The search probability of ten-qubit with respect to the iteration number under the cavity decay rate (a) 
/Lt = 0, (b) ^ = 0.05, and (c) fi = 0.1, respectively. 

Fig. 4. Infidelity of ten-qubit phase gate versus time delay in the case of r/ = 10, where fi = 0.05 and /x = 0.1 are 
plotted by solid and dashed curves, respectively. 
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